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Abstract 

We study a general class of line-soliton solutions of the Kadomtsev-Petviashvili II (KPII) equation by 
investigating the Wronskian form of its tau-function. We show that, in addition to previously known 
line-soliton solutions, this class also contains a large variety of new multi-soliton solutions, many of 
which exhibit nontrivial spatial interaction patterns. We also show that, in general, such solutions consist 
of unequal numbers of incoming and outgoing line solitons. From the asymptotic analysis of the tau- 
function, we explicitly characterize the incoming and outgoing line-solitons of this class of solutions. 
We illustrate these results by discussing several examples. 



1 Introduction 

The Kadomtsev-Petviashvili (KP) equation 

s(- 4 » + a? +6 "a;) +to V = <L1) 

where u = u(x,y,t) and a 2 = ±1, is one of the prototypical (2+l)-dimensional integrable nonlinear partial 
differential equations. The case G 2 = — 1 is known as the KPI equation, and a 2 = 1 as the KPII equa- 
tion. Originally derived [12] as a model for small-amplitude, long- wavelength, weakly two-dimensional 
(y-variation much slower than the x-variation) solitary waves in a weakly dispersive medium, the KP equa- 
tion arises in disparate physical settings including water waves and plasmas, astrophysics, cosmology, optics, 
magnetics, anisotropic two-dimensional lattices and Bose-Einstein condensation. The remarkably rich math- 
ematical structure underlying the KP equation, its integrability and large classes of exact solutions have been 
studied extensively for the past thirty years, and are documented in several monographs lTll4ll9l fT5lfT8ll20l . 

In this article we study a large class of solitary wave solutions of the KPII equation. It is well-known 
(e.g., see Refs. |3[l5l) that solutions of the KPII equation can be expressed as 

a 2 

u(x,y,t) =2—\ogx(x,y,t), (1.2) 
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where the tau function x(x,y,t) is given in terms of the Wronskian determinant 181 1131 



x(x,y,t) =Wr(f h ...,fy 
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with /W = d'f/dx', and where the functions fi,...,fy are a set of linearly independent solutions of the 
linear system 

d JL = d ll d JL = d Il 

dy dx 2 ' dt dx 3 

It should be noted that Eq. dl.3t can also be obtained as the composition of N Darboux transformations 
for KPII [ 15 ]. Note also that the Lax pair of the KP equation is given by [5 1 od y f — d 2 f + uf = and 
d t f — d^f + 6u(d x u)f + 3c(d~ l d y ii)f = 0. Thus, the functions fi,...,fy in Eqs. d 1 .3I > are precisely N 
solutions of the zero-potential Lax pair of KPII. A one-soliton solution of the KPII equation is obtained 
by choosing N = I and f(x,y,t) = e 6 ' +e® 2 , where 



Q m (x, y,t)= k m x + klj + k^t + 8 m , 



(1.5) 



with k m ,Q m fi 6 M, m = 1,2 and with k\ ^ k 2 for nontrivial solutions. Without loss of generality, one can 
order the parameters as &i < k 2 . The above choice yields the following traveling- wave solution 



u(x,y,t) = 5(^2-^1) sech 2 2(62-61) = <£(k-x + Cflf) , 
where x = (x,y). The wavevector k = (l x ,l y ) and the frequency go are given by 



k = (k\ —k 2 ,kf -kj) , 00 = ^-^2, 



(1.6) 



(1.7) 



and they satisfy the nonlinear dispersion relation 

-4(ol x + 4 + 31^ = 0. (1.8) 

The solution in Eq. dl.6t is localized along points satisfying the equation 61 =62, which defines a line in 
the the xy-plane, Such solitary wave solutions of the KPII equation are thus called line solitons. They are 
stable with respect to transverse perturbations. It is worth mentioning here that the KPI equation (namely, 
Eq. dl.lt with a 2 = — 1) also admits line-soliton solutions, but these solutions are not stable with respect to 
small tranverse perturbations. 

Equation dl.6t also implies that, apart from a constant 610 — 62,0 (corresponding to an overall translation 
of the solution), a line soliton of KPII is characterized by either the phase parameters k\,k 2 , or by two 
physical parameters, namely, the soliton amplitude a and the soliton direction c, defined respectively as 

a = k 2 — k\, c = k\+k 2 . (L9) 

Note that c = tana, where a is the angle, measured counterclockwise, between the line soliton and the 
positive y-axis. Hence, the soliton direction c can also be viewed as the "velocity" of the soliton in the 
xy-plane: c = —dx/dy = l y /l x . For any given choice of amplitude and direction of the soliton, one obtains 
the phase parameters k\ 2 uniquely as k\ = \{c — a) and k 2 = ^(c + a). Finally, note that when c = 
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(equivalently, k\ = — fe), the solution in Eq. (I1.6t becomes ^-independent and reduces to the one-soliton 
solution of the Korteweg-de Vries (KdV) equation. 

Similar to KdV, it is also possible to obtain multi-soliton solutions of the KPII equation. As y — > ±00, 
each of these multi-soliton solutions consists of a number of line solitons which are exponentially separated, 
and are sorted according to their directions, with increasing values of c from left to right as y — > —00 and 
increasing values of c from right to left as y — > 00. However, the multi-soliton solution space of the KPII 
equation turns out to be much richer than that of the (l+l)-dimensional KdV equation due to the dependence 
of the KPII solutions on the additional spatial variable y. 

It is possible to construct a general family of multi-soliton solutions via the Wronskian formalism of 
Eq. (I1.3t by choosing M phases 61 , ... ,6m defined as in Eq. dl.5t with distinct real phase parameters k\ < 
&2 < • • • < k,M and then considering the functions f\,...,fyva Eq. dl.3t defined by 

m 

f„(x,y,t) = Z V e > n = l,2,...,N, (1.10) 

m= 1 

The constant coefficients a n m define the N x M coefficient matrix A := (a nM1 ), which is required to be of 
full rank (i.e., rank(A) = AO and all of whose non-zero N x N minors must be sign definite. The full rank 
condition is necessary and sufficient for the functions /„ in Eq. dl.lOt to be linearly independent. The sign 
definiteness of the non-zero minors is sufficient to ensure that the tau function z(x,y,t) has no zeros in the 
xy-plane for all t, so that the KPII solution u(x,y,t) resulting from Eq. dl.2t is non-singular. 

One of the main results of this work (Theorem \3.6l is to show that, when the coefficient matrix A 
satisfies certain irreducibility conditions (cf. Definition l2.2t . Eq. (ll.lOt leads to a multi-soliton configuration 
which consists of N- asymptotic line solitons as y — > —00 and N + asymptotic line solitons as y — > 00, with 
N- =M—N and AO- = N, and where each of the asymptotic line solitons has the form of a plane wave similar 
to the one-soliton solution in Eq. (II ,6b . We refer to these multi-soliton configurations as the (N-,N+)- 
soliton solutions of KPII; also, we will call incoming line solitons the asymptotic line solitons as y — ► —00 
and outgoing line solitons those as y — > 00. The amplitudes, directions and even the number of incoming 
solitons are in general different from those of the outgoing ones, depending on the values of M, N, the 
phase parameters k\, . . . ,&m and the coefficient matrix A. Moreover, these multi-soliton solutions of KPII 
exhibit a variety of spatial interaction patterns which include the formation of intermediate line solitons and 
web structures in the xy-plane (2l|3]OlEIll23- In contrast, for the previously known ll23l l6l fT5l ordinary 
soliton solutions of KPII (cf. section |4} and solutions of KdV the solitons experience only a phase shift 
after collision. In several cases studied so far, the existence of these nontrivial spatial features was found 
to be related to the presence of resonant soliton interactions IH71IT91I2T1 . Several examples of these novel 
(N-,N + ) -soliton solutions of KPII are discussed throughout this work (e.g., see Figs. ^4}. 

If M = 2N, it follows that N- = N + = N, i.e., the numbers incoming and outgoing asymptotic line 
solitons are the same; we call the resulting solitons the N -soliton solutions of KPII. Among these, there is 
an important sub-class of solutions, for which the amplitudes and directions of the outgoing line solitons 
coincide with those of the incoming line solitons; we call these the elastic Af-soliton solutions of KPII. Elastic 
Af-soliton solutions possess a number of interesting features of their own, and their specific properties are 
further studied in Refs. I31 IT31 . 

We note that multi-soliton solutions exhibiting nontrivial spatial structures and interaction patterns were 
also recently found in other (2+l)-dimensional integrable equations. For example, solutions with soliton 
resonance and web structure were presented in Refs. fTOllTTl for a coupled KP system, and similar solutions 
were also found in Ref . [ 14 ] in discrete soliton systems such as the two-dimensional Toda lattice, together 
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with its fully-discrete and ultra-discrete analogues. In other words, the existence of these solutions appears 
to be a rather common feature of (2+l)-dimensional integrable systems. Thus, we expect that the scope of 
the results described in this work will not be limited to the KP equation alone, but will also be applicable to 
a variety of other (2+l)-dimensional integrable systems. 

2 The taii-function and the asymptotic line solitons 

In this section we investigate the properties of the tau-f unction in Eq. dl.3l > when the N functions fi,.-.,fy 
are are chosen according to Eq. dl.lOt as linear combinations of M exponentials e 9 ' , . . . ,e 9M . We should 
emphasize that Eq. (11.101) represents the most general form for the functions involving linear combinations of 
exponential phases. Since the elements of the N xM coefficient matrix A = (a n ^ m ) are the linear combination 
coefficients of the functions f\,...,fn, one can naturally identify each /„ with one of the rows of A and each 
phase Q m with one of the columns of A, and viceversa. In this section we examine the asymptotic behavior 
of the tau-function in the xy-plane as y — » ±00. It is clear that, with the above choice of functions, the 
tau-function is a linear combination of exponentials. Consequently, the leading order behavior of the tau- 
function as y — > ±00 in a given asymptotic sector of the xy-plane is governed by those exponential terms 
which are dominant in that sector. A systematic analysis of the dominant exponential phases allows us to 
characterize the incoming and outgoing line solitons of (A^_,A^ + )-soliton solutions of KPII. 

2.1 Basic properties of the tau-function 

We start by presenting some general properties of the tau-function. Without loss of generality, throughout 
this work we choose the phase parameters k m to be distinct and well-ordered as k\ < £2 < • • • < 

Lemma 2.1 Suppose %n,m = Wr(/j , . . . as in Eq. (U.3t . with the functions f\ , . . . , fy given by Eq. dl.lOt . 
Then 

z NM (x,y,t)=det(A®K T ), (2.1) 

where A = (a n , m ) is the N xM coefficient matrix, & = diag(e 01 , . . . ,e® M ), and the N xM matrix K is given 
by 

/ 1 1 •■• 1 \ 

K= . . . , 

where the superscript T denotes matrix transpose. Moreover, %n,m can be expressed as 

?N,M(x,y,t) = E V(m 1 ,...,m N )A(mi,...,m N ) exp[ 6 OTWBAr ] , (2.2) 

l<mi <mi<---<n\N <M 

where Q mit ...m N denotes the phase combination 

®m l ,...,m N (x,y,t) = Q m (x,y,t) -\ ^Q mN (x,y,t) , (2.3) 

A(m\ ,. .. ,/?Jjv) denotes the N x N minor of A obtained by selecting columns mi,. .. ,m^, and V(m\ ,. .. ,m^) 
denotes the Van der Monde determinant 

V(m 1 ,...,mi f )= Yl ( km n ~ km n ) ' ^ 2A ^ 

l<s { <s 2 <N 
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Proof. Equation (I2.lt follows by direct computation of the Wronskian determinant d 1 .3I >- Next, to prove 
Eq. (12.21) apply the Binet-Cauchy theorem to expand the determinant in Eq. (I2.lt and note that the N x 
N minor of K obtained by selecting columns 1 <m\ < ■■■ < < M is given by the Van der Monde 
determinant V(m\ m^). □ 

From Lemma lXTl we have the following basic properties of the tau-f unction: 

(i) The spatio-temporal dependence of the tau-function in Eq. (I2.2t is confined to a sum of exponential 
phase combinations 8 mi ,..., mjv which according to Eq. (I2.3t are linear in x, y,t. Moreover, all the Van der 
Monde determinants V(m\, . . . ,m^) are positive, as the phase parameters ki, . . . ,1cm are well-ordered. 
Note that a sufficient condition for the tau-function (I2.2t to generate a non-singular solution of KPII 
is that it is sign-definite for all (x,y,t) 6K 3 . In turn, a sufficient condition for the tau-function (12.21) 
to be sign-definite is that the minors of the coefficient matrix A are either all non-negative or all 
non-positive. Note however that it is not clear at present whether these conditions are also necessary. 

(ii) Each exponential term in the tau-function of Eq. (12.21) contains combinations of N distinct phases 
8 m , , . . . ,Q,„ N identified by integers mi,... chosen from {1, .. . ,M}. Thus, the maximum num- 
ber of terms in the tau-function is given by the binomial coefficient (^) . However, a given phase 
combination Q mij ,.. >mN is actually present in the tau-function if and only if the corresponding minor 
A(mi,. . . ,m#) is non-zero. 

(iii) If M < N the functions /i, ...,/# are linearly dependent; in this case there are no terms in the summa- 
tion in Eq. (12.2b . and therefore the tau-function %N,M{x,y,t) is identically zero. Also, if M = N, there 
is only one term in the summation (corresponding to the determinant of A); then %Njt(x,y,t) depends 
linearly on x and therefore it generates the trivial solution of KP. Finally, if rank(A) < N, all N xN 
minors of A vanish identically, leading to the trivial solution tffjf(x,y,t) = 0. Therefore, for nontrivial 
solutions one needs M > N and rank(A) = N. 

(iv) The transformation A — > A' = GA with G € GL(iV,M) (corresponding to elementary row operations 
on A) amounts to an overall rescaling x(x,y,t) — > x'(x,y,t) = det(G)z(x,y,t) of the tau-function (I2.lt . 
Such rescaling leaves the solution u{x,y,t) in Eq. (I1.2t invariant. This reflects the fact that N indepen- 
dent linear combinations of the functions fi,...,fy in Eq. dl.lOt generate equivalent tau-f unctions. 
This GL(iV,R) gauge freedom can be exploited to choose the coefficient matrix A in Eq. (I2.lt to be 
in reduced row-echelon form (RREF). As is well-known, the GL(Af,R) invariance means that the 
tau-function (I2.lt represents a point in the real Grassmannian Gr(iV,M). 

(v) Suppose that one of the functions in Eq. (ll.lOt contains only one exponential term; that is, suppose 
f p = a Piq e Qq with a pjn = Mm ^ q. Then it is A(m\, . . . ,m^) = whenever q ^ {mi, . . . ,otjv}, and the 
resulting tau-function (I2.2t can be expressed as lNju(x,y,t) = e B ix'(x,y,t), where i'(x,y,t) is a linear 
combination of exponential terms containing combinations of N — I distinct phases chosen from the 
remaining M —I phases (that is, all M phases but 6 ¥ ). From Eq. (I1.2t it is evident that Tffjn(x,y,t) and 
z'(x,y,t) generate the same solution of KP. Moreover, the function x'(x,y,t) is effectively equivalent 
to a tau-function Tn-im-i (x,y,t) with a coefficient matrix obtained by deleting the p-th row and q-th 

column of A. Hence, the tau-function Tffjt(x,y,t) is reducible to another tau-function x^v— i.a/ l (x,y,t) 

obtained from a Wronskian of N — I functions with M — 1 distinct phases. 

In accordance with the above remarks, throughout this work we consider the coefficient matrix A to be 
in RREF. Also, to avoid trivial and singular cases, from now on we assume that M > N and rank(A) = N, 
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and that all non-zero N x N minors of A are positive. Finally, we assume that A satisfies the following 
irreducibility conditions: 

Definition 2.2 (Irreducibility) A matrix A of rank N is said to be irreducible if, in RREF: 

( i) Each column of A contains at least one non-zero element. 

( ii) Each row of A contains at least one non-zero element in addition to the pivot. 

Condition (i) in Definition 12. 21 requires that each exponential phase appear in at least one of the functions 
f\,...,fy; condition (ii) requires that each function contains at least two exponential phases. The reason for 
condition (i) should be obvious, for if A contains a zero column, the corresponding phase is absent from the 
tau-function, which can then be re-expressed in terms of an irreducible N x (M — 1 ) matrix. The reason for 
condition (ii) is to avoid reducible situations like those in part (v) of the above remarks. Note also that if an 
N x M matrix A is irreducible, then M > N. 

2.2 Dominant phase combinations and index pairs 

We now study the asymptotic behavior of the tau-function in the xy-plane for large values of \y\ and finite 
values of t. Let denote the set of all phase combinations 6 mi ,..., m/v such that A(m\, . . . ^ 0, that is, the 
set of phase combinations that are actually present in the tau-function x(x,y,t). 

Definition 2.3 (Dominant phase) A given phase combination Q mi mN E is said to be dominant for the 
tau-function x(x,y,t) ofEq. i2.2t in a regionR E M 3 ifQ m > lr ,. jm > N (x,y,t) < B mu ..., mN (x,y,t) for all Q m ' v .,. ;m ' E & 
and for all (x,y,t) E R. The region R is called the dominant region ofB mii ,_, mN . 

As the phase combinations B mi ,...,m N (x,y,t) are linear functions of x,y and t, each of the inequalities in 
Definition 12 . 31 defines a convex subset of R 3 . The dominant region R associated to each phase combination 
is also convex, since it is given by the intersection of finitely many convex subsets. Furthermore, since 
the phase combinations are defined globally on M 3 , each point (x,y,t) E M 3 belongs to some dominant 
region 7?. As a result, we obtain a partition of the entire M 3 into a finite number of convex dominant regions, 
intersecting only at points on the boundaries of each region. It is important to note that such boundaries 
always exist whenever there is more than one phase combination in the tau-function, because then there are 
more than one dominant region in IR 3 . The significance of the dominant regions lies in the following: 

Lemma 2.4 The solution u(x,y,t) of the KP equation generated by the tau-function (11.31) is exponentially 
small at all points in the interior of any dominant region. Thus, the solution is localized only at the bound- 
aries of the dominant regions, where a balance exists between two or more dominant phase combinations in 
the tau-function ofEq. \2.2\ . 

Proof. Let R be the dominant region associated to 8 m , OTJ¥ , which is therefore the only dominant phase in 
the interior of R. Then from Eq. \2.2\ we have that % .M(x,y, t) ~ 0(e 6m i>~< m N ) in the interior of R. As a result, 
log TNjf(x,y,t) locally becomes a linear function of x apart from exponentially small terms. Then it follows 
from Eq. (11.21) that the solution u(x,y,t) of KP will be exponentially small at all such interior points. □ 

The boundary between any two adjacent dominant regions is the set of points across which a transition 
from one dominant phase combination 8 m|v .. jmjv to another dominant phase combination Q m i m ' N takes 
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place. Such boundary is therefore identified by the equation 6 m , „, N = 6 m / r .. jOT ' , which defines a line in the 

xy-plane for fixed values of t. The simplest instance of a transition between dominant phase combinations 
arises for the one-soliton solution ( U.6t . which is localized along the line 61 = 62 defining the boundary 
of the two regions of the xy-plane where 61 and 62 dominate. In the one-soliton case, these two regions 
are simply half -planes, but in the general case the dominant regions are more complicated, although the 
solution u(x,y,t) is still localized along the boundaries of these regions, corresponding to similar phase 
transitions. For example, Fig. ^ illustrates a (2, l)-soliton known as a Miles resonance fFTll (also called 
a Y-junction), generated by the tau-function T12 = e 6 ' +e 02 +e 93 . In this case, the xy-plane is partitioned 
into three dominant regions corresponding to each of the dominant phases 61, 62 and 63. Once again, the 
solution u{x,y,t) is exponentially small in the interior of each dominant regions, and is localized along the 
phase transition boundaries: here, 61 = 62, 61 = 63 and 62 = 63. It should also be noted that some of these 
regions have infinite extension in the xy-plane, while others are bounded, as in the case of resonant soliton 
solutions, described in section|4]and Ref. [2 1. Each phase transition which occurs asymptotically as y — > ±00 
defines an asymptotic line soliton, which is infinitely extended in the xy-plane. 

When studying the asymptotics of the tau-function for large |y| it is useful to employ coordinate frames 
parametrized by the values of direction c. That is, we consider the limit y — > ±00 along the straight lines 

L c :x + cy = ^. (2.5) 

Note that c increases counterclockwise, namely from the positive x-axis to the negative x-axis for y > and 
from the negative x-axis to the positive x-axis for y < 0. From Eqs. dl.5t and (I2.5t . the exponential phases 
along L c are 6 m = k m (k m — c) y + k m £, + kf n t + 6„,.o. The difference between two such phases along L c is then 
given by 

9 m - 6m' = {k„ - k m < ) {k m + k m < - c)y + {k m - k m ' %+{k m - k m f)t + Q m ,0 - 8m',0 , (2.6a) 
and the difference between any two phase combinations along L c is given by 

Qm u ...,m N -9m', ,...,m' N = ^ E (kmj ~ Ki') {k mj + krf. - c) ^ y + 8(£,f) , (2.6b) 

where 8(^,?) = lI=i \{kmi ~Ki')^ + — k 3 , )t + 9 m ,.o — 9m'. ol ■ m particular, the single-phase-transition 

J L J j J J ' 

line L m m i : 6 m = Q m i, which will play an important role below, is given by Eq. (I2.5t with c m ,„/ = k m + k m >. 

Before proceeding further, we introduce the following notations which will be employed throughout this 
article. We denote by A [m] £ M. N the m-th column of A, and we denote by A [m\ , ... , m r ] the N xr submatrix 
obtained by selecting the r columns A [mi], . . . ,A[m r ]. We also label the ,/V pivot columns of an irreducible 
N xM coefficient matrix A by A[e{\, . . . ,A[<?jv], with 1 = e\ < <?2 < ■ • • < < M, and we label the M — N 
non-pivot columns by A[gi], . . . ,A[gM-/v], where 1 < g\ < g2 < ■ ■ ■ < gM-N = Af. Note that A has pivot 
columns because it is rank A^; also, e\ = \ since A is in RREF, and e^ < M since it is irreducible. We now 
establish a result that will be useful in order to characterize the asymptotics of the tau-function. 

Theorem 2.5 (Single-phase transition) Asymptotically as y — > ±°°, and for generic values of the phase 
parameters k\ , . . . , the dominant phase combinations in the tau-function (I2.2t exhibit the following be- 
havior in the xy-plane: 

(i) the set of dominant phase combinations remains invariant in time for finite values oft. 
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( ii) the dominant phase combinations in any two adjacent dominant regions contain N — 1 common 
phases. 

The proof of Theorem 12. 5 l is given in the Appendix. 

Consider the single-phase transition as y — > ±oo in which a phase 6, from the dominant phase combina- 
tion in one region is replaced by another phase Qj to produce the dominant phase combination in the adjacent 
region. We refer to this transition as an i —> j transition, which takes place along the line Lu : 6,- = Qj whose 
direction in the xy-plane is given by ctj = kj + kj. As y — > °°, it is clear from Eq. (I2.6at that, if k{ < kj, the 
transition i — > j takes place from the left of the line L,j to its right, while if kj > kj the transition i — » j takes 
place from the right of the line L, <j to its left. Thus, as y — » °°, each dominant phase region /? is bounded 
on the left by the transition line Ly given by to the minimum value of cy that corresponds to an allowed 
transition, and, similarly, on the right by the transition line Ly given by the maximum value of c,-j that 
corresponds to an allowed transition. Here, an allowed transition from one dominant phase combination to 
another means that the minors associated with those phase combinations in the tau-function of Eq. J2.2t . 
are both non-zero. In turn, these non-vanishing minors determine the values of cu corresponding to the 
allowed single-phase transitions. A similar statement can be made for transitions occurring as y —> — oo. So, 
each dominant phase region R as y — > ±°o has boundaries defined by a counterclockwise and a clockwise 
single-phase transitions which can be identified as follows: 

Corollary 2.6 Suppose that 6 mij ..., mA , is the dominant phase combination on a region R asymptotically as 
y — > ±oo. Let J be the complement of the index set {m\,m2, ■ ■ ■ ,m^} in {1,2, .. . ,M}. Also, for each element 
j € J, define a corresponding L C {mi ,rri2,---, m^} as the set of all indices m r € {m\ ,rri2, ■ ■ ■ , m^} such that 
the minor A(m\, . . . ,m r -i,j,m r +i,. . . ,itin) 0- Then: 

(i) as y — > oo, the directions of the counterclockwise and clockwise transition boundaries ofR are respec- 
tively given by 

c + = min [cij] with > kj, c_ = max [a j] with kj <kj. (2.7a) 

ieljJeJ ' ieljJeJ 

( ii) as y — > — °°, the directions of the counterclockwise and clockwise transition boundaries of R are 
respectively given by 

c + = min [en] with kj<kj, c_ = max [a j] with (2.7b) 

ieijjeJ ' ' ieljJeJ 

The results of Theorem 12.51 and Corollary I2.6l can now be used to determine the asymptotic behavior of 
the tau-function of Eq. dl.3t . thereby obtaining an important characterization of the asymptotic line solitons 
corresponding to (Af_,/V + )-soliton solutions of the KPII equation. Namely, for the tau-function lN^t(x,y,t) 
of Eq. (I2.2t with generic values of the phase parameters k\, . . . ,kw we have the following: 

(i) As y —> ±oo, the dominant phase combinations of the tau-function in adjacent regions of the Ay-plane 
contain N — I common phases and differ by only a single phase. The transition between any two 
such dominant phase combinations 6i, m2r .. !mjv and Qj. m2r .. MN occurs along the line L,j : 6 ; - = 8 7 -, where 
a single phase 8, in the dominant phase combination is replaced by a phase Qj. Moreover, if the 
dominant phase combination Qi. mi ,...,m N in a given region is known, the transition line L, y - and the 
dominant phase combination 0y, m2 ,..., mjv are determined via Corollary 12.61 In particular, Eqs. (I2.7t 
for c± determine explicitly the pair of phase parameters k[ and kj corresponding to the single-phase 
transition i — ► j across each boundary Lij of a given dominant phase region. 
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(ii) As y — ► ±00 along the line L,-j, the asymptotic behavior of the tau-function is determined by the 
balance between the two dominant phase combinations Qi : m 2 ,...,m N and Qj : m 2 ,...,m N , and is given by 

x NM (x,y,t) ~V(i,m 2 , . . .,m N )A(i,m 2 , . . .,m N )e Q '^ «w +V(j,m 2 , . . .,m N )A(j,m 2 , . . .,m N )e Q ^ »w , 

(2.8a) 

where V(mi,. . . ,mjy) is the Van der Monde determinant defined in Eq. (12.41) . and where the minors 
A(i,m 2 ,...,niN) and A (j,m 2 , ... ,mjv) of the coefficient matrix A are both non-zero. The solution 
«(x,;y,f) of the KPII equation in a neighborhood of such a single-phase transition is then obtained 
from Eq. dl.2t as, 

u(x,y,t) ~ -* ; -) 2 sech 2 [±(9; - 6;)] • (2.8b) 

Moreover, Lemma I2~4l and Theorem 12.51 together imply that the solution of the KPII equation is ex- 
ponentially small everywhere in the xy-plane except at the locations of such single-phase transitions. 
Equation d2.8bt . which is a traveling wave solution satisfying the dispersion relation in Eq. dl.8t . co- 
incides with the one-soliton solution in Eq. dl.6t . Thus, it defines an asymptotic line soliton associated 
with the single-phase transition i — ► j. The phase parameters k { and kj associated with the single-phase 
transition i — > j are determined by Eqs. (12.71) : the soliton amplitude is thus given by a;j = \kj —kj\, 
and the soliton direction is given by the direction of L, ; , which is cu = ki + kj. 

(iii) All of the asymptotic line solitons resulting from single-phase transitions such as the one described 
above are invariant in time, in the sense that their number, amplitudes and directions are constants. 

Motivated by these results, we label each asymptotic line soliton by the index pair [i,j] which uniquely 
identifies the phase parameters ki and kj in the ordered set {k\, . . . ,&m}- The results summarized in the 
above remarks can be applied to explicitly delineate the dominant phase combinations and the asymptotic 
line solitons associated with the tau-function of a given (N-,N + ) -soliton solution of the KPII equation, as 
illustrated by the following example. 

Example 2.7 When N = 2 and M = 4, Lemma l2~T1 implies that the tau-function x(x,y,t) is given by 

z(x,y,t)=Wr(f l ,f 2 )= I {k ml -k m )A{m,m')^- + ^\ (2.9) 

l<m<m'<4 

where the four phases are given by 6„, = k m x + k^y + k^t + 6,„ for m = 1 , . . . , 4, as in Eq. (11.51) . and where 
the phase parameters are ordered as fci < ■ • • < £4. We consider the line-soliton solution constructed from 
the two functions f\ = e 6 ' + e 02 and f 2 = e 03 + e 94 , so that the associated 2x4 coefficient matrix is 

a A 1 (A 

A= (o 1 ij- (2 - 10) 

Then A (1,2) = A (3, 4) =0, and the remaining four minors are all equal to one. We apply Corollary 12.61 
to determine the asymptotic line solitons associated with the tau-function in Eq. i2.9i . First note from the 
expression Q m m i = (k m + k m i)x + kf n + k 2 m i)y + (k m + kf n ,)t + m o + Q m ',0 that for every finite value of y the 
dominant phase combination as x — » —00 is given by 81 3, which corresponds to the minimum value of k m + 
k m i such that A(m, m') 7^ (cf. Definition 12. 3 l l. Let us denote by #13 the region of the xy-plane where 6^3 is 
the dominant phase. The transition boundaries of R\ 3 are determined by applying Corollary I2.6l as follows: 
The complement of the index set {1,3} is / = {2,4}. When j = 2 G J, we have A(l,2) = but A(2,3) 7^ 0; 
hence I 2 = {1}. Similarly, when j = 4 we have I 4 = {3} because A (1,4) / but A(4,3) = 0. Thus the 
possible transitions i — > j from R\ 3 are 1 — > 2 and 3 — > 4. As y — > 00, the second of Eqs. i2.lt implies that 
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Figure 1 : Dominant phase combinations in the different regions of the xy-plane (labeled by the indices in parentheses) 
and the asymptotic line solitons (labeled by the indices in square braces) for two different line soliton solutions: (a) a 
fundamental Miles resonance (Y-junction) produced by the tau-function with N = l,M = 3and (&i,&2,&3) = (— 1,0, 5) 
at t = 0; (b) an ordinary two-soliton solution, produced by the coefficient matrix in Example 12 . 71 with (k\,... ,£4) = 
(— |, — i,0, 1) at t = (see text for details). Here and in all of the following figures, the horizontal and vertical axes 
are respectively x andy, and the graphs show contour lines of log u(x,y,t) at a fixed value of t. 



the clockwise transition boundary of /?i 3 is given by the transition line L3.4, whose direction C3 4 = £3 + ^4 
is greater than the direction c\ 2 = k\ + k 2 of the line L\ 2 . Across the transition line L 3 4 , the dominant phase 
combination switches from 61 3 to 614, onto the corresponding dominant region, which we denote R\a- 
Similarly, as y — ► —00, the first of Eqs. (I2.7bt implies that the counterclockwise transition boundary of 
is given by the transition line L\ 2 , whose direction is less than the direction C3.4 of the line L3.4. This 
implies that the dominant phase combination and dominant region change to 62.3 and /?2.3, respectively. 
Applying Corollary 12. 6l again to the region R 2 .z as y — > — °°, one finds J = {1,4} with l\ = {2} and I4 = {3}, 
so the possible transitions from R 2 3 are 2 — > 1 and 3^4. The 2 — ► 1 transition corresponds to a clockwise 
transition from R 2 i, back to R\ i, whereas the 3^4 transition corresponds to a counterclockwise transition 
from R 2 t, to the region R 2 ,4, where 62,4 is the dominant phase combination. Continuing counterclockwise 
from Rip we finally obtain the following dominant phase regions asymptotically as y — > ±00, together with 
the associated single-phase transitions: 

1— >2 3^4 2— >1 4^3 
^1,3 ^2,3 ^3,4 >Rl,A ^1,3 • (2.11) 

It is then clear that there are two asymptotic line solitons as y — > —00 as well as y — > 00, and in both cases 
they correspond to the lines 61 = 62 and 63 = 64. The dominant phase regions, denoted by indices (m,m'), 
and the asymptotic line solitons, identified by the index pairs [i,j], are illustrated in Fig.^>. 

In the following section we obtain several results that will allow us to identify more precisely the index 
pairs corresponding to each asymptotic line soliton. In addition, we will prove a general result concerning the 
numbers of asymptotic line solitons present in any (N-,N + ) -soli ton solution corresponding to a tau-function 
with an arbitrary number of functions fi,...,fy and arbitrary linear combinations of the exponential phases 
e Q[ , . . . ,e® M in each function. 
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3 Asymptotic line solitons and the coefficient matrix 



In this section we continue our investigation of the tau-function in the general setting introduced in sectional 
We have seen in the previous section that an asymptotic line soliton corresponds to a dominant balance 
between two phase combinations in the tau-function. But we still need to identify which phase combinations 
in a given tau-function are indeed dominant as y — > ±00. This requires a detailed study of the structure of the 
N xM coefficient matrix A associated with the tau-function. In this section we carry out this analysis, which 
enables us to explicitly identify all the asymptotic line solitons of a given tau-function in an algorithmic 
fashion. One of our main results of this section will be to establish that, for arbitrary values of N and M, and 
for irreducible coefficient matrices (cf. Defmition l2.2t with non-negative N xN minors, the tau function dl .3I > 
produces an (N-,N + ) -soli ton solution with = M — N and N + = N, i.e., a solution in which there are 
N- = M — N asymptotic line solitons as y — > —00 and N + = N asymptotic line solitons as y — > 00. 

3.1 Dominant phases and structure of the coefficient matrix 

We begin by presenting a simple yet useful result that will be frequently used to determine the dominant 
phase combinations in the tau-function as y —> ±00. 

Lemma 3.1 (Dominant phase conditions) As y — > ±00 along the line I4 j : 9, = dj with i < j, the exponential 
phases 61 , . . . , Qm satisfy the following relations. 

(i) As y — > 00, Q m < 0^ Vm G {i+ 1, ... , j — 1}, and m > 6*, Vm G {1, . .. , i — l,j + 1, .. . ,M}, where 
0* := 0; = 0y. 

(ii) As y — > -00, 6 m > 6*, Vm G {i+ 1, . .. , j — 1}, while m < 0*, Vm G {1, — 1J+ 1, .. . ,M}. 



Proof. It follows from Eq. (I2.6at that, along the line L,j, the difference beetween any two exponential 
phases 0,„ and 0,,,/ is given by 

0m -0,„' = {k m -k ml ){(k m +k ml )-{ki + kj)\y + h'{^t), (3.1) 

where 8'(£,f) is a linear function of £, and t and which also depends on the constants m 0, m ' o, 0i.o and 0;o, 
and where we used the fact that the direction of the line Ljj is c ; j = fej +kj. It is clear that the sign of 
0m — m ' as y — > ±00 and for finite values of ^ and t is determined by the coefficient of y in the right-hand 
side of Eq. (I3.lt Then, setting m' = i (or m' = j) in Eq. (I3.lt one obtains the desired inequalities. □ 

Lemma l3~Tl which is illustrated in Fig.|2j will be used to obtain a set of conditions that are necessary 
for a given pair of phase combinations in the tau-function to be dominant. These conditions are given in 
terms of the vanishing of certain N x N minors of the coefficient matrix A, and they determine which phase 
combinations are present (or absent) in the tau-function of Eq. dl -3b - In order to derive these conditions, it is 
convenient to introduce two submatrices Pjj and Qij associated with any index pair [/, j] with 1 < i < j < M, 
and given by 

/'// /\ 1.2 / I../ - 1 M . Q Lj M ■ \ ./ 1 . (3.2) 

The matrix P, -j contains the consecutive columns of A to the left of column A [i] and those to the right of 
column A[j], while Qij contains the consecutive columns of A between columns A[i] and A[j}. Using the 
matrices and Q/j and the dominant phase conditions in Lemma 13171 we then have: 
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As y — oo along 6, = 9- 



kf k 2 k,_i k, k j+ i kj_^ kj kj +1 ■ ■ ■ /c M _ 1 k M k 



e m >e / =e 7 Wm<w ' e m >e / =e y 

As y —» - oo along 0, = 



ky k 2 kj_-\ kj /c /+1 kj_-\ k } k M _ A k M k 

\ j 

\ J V V J 

Y Y 

e m <e,= e y e m > e /= 9 / e m <e,= e y 

Figure 2: Relations among the exponential phases as y — > ±°° along the direction L,j : 8, = 9y. 



Lemma 3.2 (Vanishing minor conditions) Suppose that the index pair [i,j] identifies an asymptotic line 
soliton. Let the two dominant phase combinations along the line Ljj : 6 ( = 6 ; be given by 

Qi, PU ~,Pr,qi,-,<is and ®j,Pi,-,Pr,qi,...,q s > and let Mh Pi > ■ ■ ■ i Pr , <?i , • • • , , p x , . . . , p r , q\ , . . . , q s ) be the 
corresponding non-zero minors where A[p\], . . . ,A[p r ] € Pjj and A[q{\,. . . ,A[q s ] E Qij. 

(i) If identifies an asymptotic line soliton as y — > oo, then 

(a) allN xN minors obtained by replacing one of the columns A[i],A[j],A[qi], . . . ,A[q s ]from either 
A(i,pi,...,p r ,q\,... ,q s ) or A(j,pi,. ..,p r ,q\,...,q s ) with any column A[p] € are zero; 

(b) all N xN minors obtained by replacing one of the columns A[qi], . . . ,A[q s ] from either 
A(i,pi, . . . , p r ,qi , ■ ■ ■ ,q s ) or A(j,pi,. . . ,p r ,q\, ■ ■ ■ ,q s ) with either A[i] orA[j], are zero. 

(ii) If [i,j] identifies an asymptotic line soliton as y — > — oo, then 

(a) allNxN minors obtained by replacing one of the columns A i .A j .A p\ .. . . ,A[p r ]from either 
A(i,pi,...,p r ,qi,... ,q s ) or A(j,p\,. . . ,p r ,qi, ■ ■ ■ ,q s ) with any column A[q] £ Qj j, are zero; 

(b) all N x N minors obtained by replacing one of the columns A[pi], . . . ,A[p r ] from either 
A(i,p\, . . . ,p r ,q\,. ..,q s ) or A(j,pi,. . . ,p r ,q\,. . . ,q s ) with either A[i] orA[j], are zero. 

Proof. All of the above conditions follow from the repeated use of the dominant phase conditions in 
Lemma l3~T1 For example, as y — > oo along the line L, j, Lemma l3~T1 implies Q p > 6 m for all p G { 1 — 
1, j + 1, . . . ,M} and for all m € {i,j,qi,...,q s }. Consequently, if condition (b) in part (i) of the Lemma 

does not hold, each of the phase combinations obtained by replacing 8 m with Q p in either Qi lPlr ., lPr> q i qs 

or ®j,p\,...,Pr,q\,-,qi wm b e g reater than both Qi lPlr .. lPrt q lr .. t q s and Qj tPlr ., tPrj q ly ., iqs . But this contradicts the 
hypothesis that &i, Pi} ... :Pr .q u ... Xh and 0y,p li ...,p, 4li ..., 9l are the dominant phase combinations as y — > oo along 
Lij. The other conditions follow in a similar fashion. □ 

We should emphasize that [i,j] denotes an asymptotic line soliton either as y —> oo or as y — > — oo. In 
general, the asymptotic solitons (and therefore the index pairs) as y — > oo and those as y — > — oo are different. 
Thus, in principle there is no relation among the matrices Pij and Qij relative to solitons as y — ► oo and those 
associated to solitons as y — > — oo. 

Lemma l3~2l allows us to determine the ranks of the submatrices Pjj and Qij associated with each asymp- 
totic line soliton This information will be exploited later in Theorem 13.61 to identify explicitly the 
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asymptotic line solitons produced by any given tau-function. The next two results are direct consequences 
of the conditions specified in Lemma l3~2l 

Lemma 3.3 (Span) Let A[p [],... ,A[p r ] be the columns from P,-; and A[qi], ... ,A[q s ] be the columns from 
Qij in the minors associated with the dominant pair of phase combinations, as in Lemma U^\ 

(i) If [i,j] identifies an asymptotic line soliton as y — > °°, the columns A[p\], . . . ,A[p r ]form a basis for the 
column space of the matrix Pj ;. 

(ii) If [i,j] identifies an asymptotic line soliton as y — > — oo, the columns A[q\], . . . ,A[q s ]form a basis for 
the column space of the matrix Qij. 

Proof. We prove part (i). Since A(i. pi y r .q\ q s ) ^ bv Lemma 13.21 the set of columns SA = 

{A[i\,A\pi],...,A\p r ],A[qi],...,A[q s ]} is a basis of R^. Hence the set {A[pi],...,A[p r ]} C SA is linearly 
independent. Moreover, for any A[p] G pj we can expand A [p] with respect to S\ : 

r s 

A[p]=aA\i}+ £ b m A[p m ] + £ c m A[q m ). (3.3) 

m= 1 m= 1 

Replacing one of the columns A [i], A [q\], . . . ,A[q s ] \nA(i,p\, . ..,p r ,q\,...,q s ) with A [p] G Py, we have from 
Lemma l3~2l i.a that 

A(p,p 1 ,...,p r ,q 1 ,...,q s ) = 0, A(i,pi,...,p r ,qi,...,q m _i,p,q m+ i,...,q s ) = 0. 

Hence in Eq. (I3.3t we have a = and c m = Vm = 1, ... ,s. Therefore A(p) G span(A[^i, . . . ,p r ]) for all 
A\p] G Pi j. Similarly, part (ii) follows from the conditions in Lemma l3~2l ii.a. □ 

Lemma 3.4 (Rank conditions) Let r be the number of columns from P, y and let s be the the number of 
columns from Qij in the minors associated with the dominant pair of phase combinations, as in Lemma U^\ 

( i) If [i, j] identifies an asymptotic line soliton asy^oo, then rank(py) = r <N — 1 and rank(P,-j |A [i] ) = 
ranki/V/A./;.) rank!/',,- A/./ } r ■ 1. 

( ii) If [i, j] identifies an asymptotic line soliton asy^ —oo, then rank((2,-j) = s < N — 1 and rank(<2,-j |A [i] ) 
= rar±(Qij\A[j])=var±(Qij\A[i,j})=s+L 

Above and hereafter, (A\B) denotes the matrix A augmented by the matrix B. 

Proof. Let us prove part (i). Since the columns A[pi], . . . ,A\p r ] form a basis for the column space of 
Pu, from Lemma l3~3l i we immediately have rank(pj) = r. Moreover, since SI = [A\i],A[p\], . . . ,A{p r ], 
A[qi], . . . ,A[q s ]} is a basis for M. N , the vectors A[i],A[pi],. . . 7 A[p r ] are linearly independent, and therefore 
rank(P,- j |A [i] ) = r + 1 . Similarly, replacing A [i] with A [j] in the previous statement we have rank(p j \A [j] ) = 
r + 1. It remains to prove that rank(P;j|A[£,y]) = r+ 1. Expanding the j-th column of A in terms of SI as in 
Lemma l3~3l we have 

A[j] = aA\i] + £ b m A[ Pm ] + £ c m A[q m ] ■ (3.4) 

m= 1 m= 1 

By replacing one of the columns A[^i], . . .A[q s ] in A(i,pi, . ..,p r ,q\,...,q s ) with A[j], from Lemma l3~2l i.b 
we have that A(i,pi, . ..,p r ,q\,.. . ,q m -i,j,qm+i, ■ ■ ■ ,qs) =0. Therefore c m = for all m = 1, . . . , s. Conse- 
quently we have A[j] G span(A[/],A[/?i], . . . ,A[/? r ]), which implies that rank(pj A[/, j]) = r + 1. Similarly, 
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using Lemma EOl ii.b one can establish the corresponding results in part (ii) for the asymptotic line solitons 
as y — > — oo. □ 



It is important to note that, even though Lemmas I3.3H3.4I were proved by using the vanishing minor 
conditions in Lemma l3~2l they provide additional information on the structure of the coefficient matrix A. 
For example, when r < N — 1 for an asymptotic line soliton as y — > oo, Lemma l3~4l vields Tanik.(Pu\A[i,j]) < 
N, and when s < N — 1 for an asymptotic line soliton as y — > — oo, Lemma l3~4l yields rank(<2,j|A[z, j]) < N. 
As a consequence, we immediately have the following additional vanishing minor conditions: 

(i) If [i, j] identifies an asymptotic line soliton as y — > oo, then 

A{i,j,p\,...,p n mi,...,m N - r -2) =0 V{mi,...,m N - r - 2 } C {1,...,M}. (3.5a) 



(ii) If identifies an asymptotic line soliton as y — > — oo, then 

A(i,j,qi,...,q s ,mi,...,m N - s - 2 ) =0 V{mi, . . . ,m N - s - 2 } C {!,..., M}. (3.5b) 



It should also be noted that, when [i,j] identifies an asymptotic line soliton as y — > oo, Lemma l3~4l i only 
provides information on Pij, and the only condition on Qij is that rank(<2ij) > s. Similarly, when 
identifies an asymptotic line soliton as y — > — oo, all we know about Pij is that rank(P/j) > r. 



3.2 Characterization of the asymptotic line solitons from the coefficient matrix 

In section 13.11 we derived several conditions that an index pair [/, j] must satisfy in order to identify an 
asymptotic line soliton. In this section we apply the results developed in section l3~Tl to obtain a complete 
characterization of the incoming and outgoing asymptotic line solitons of a generic line-soliton solution of 
the KPII equation. 



Lemma 3.5 (Pivots and non-pivots) Consider an indt 

( i) If [i, j] identifies an asymptotic line soliton asy — 
matrix A. That is, A[i] =A[e n ] with 1 <n <N. 

(ii) If identifies an asymptotic line soliton as y 
coefficient matrix A. That is, A[j] =A[g„] with 1 



x pair [i, j] with 1 < i < j <M. 

■ oo, the index i labels a pivot column of the coefficient 

— > — oo, the index j labels a non-pivot column of the 
< n<M —N. 



Proof. We first prove part (i). Suppose that 6 ! m2i ..., mA , is one of the dominant phase combinations corre- 
sponding to the asymptotic line soliton [i,j] as y — > oo. The corresponding minor A(i,m 2 , . . . ,m^) is non- 
zero. Since A is in RREF, we have A [i] = £"=l c rM e r] f° r some n <N, where e\ < ■ ■ ■ < e n < i. Therefore 
A(i,m 2 , . . . ,«iat) = £" =1 c r A(e r ,m 2 , . . . ,m^). If e„ < i, we have A[e\], . . . ,A[e n ] € where P,j is the sub- 
matrix of A defined in Eq. (13.21 . Then from condition (a) in Lemma l3~2l i we have A(e r ,m 2 , . . . ,m^) = 
Vr = 1, . . . ,n, implying that A(i,m 2 , . . . ,m#) = 0. But this is impossible, since 6;. m2i ...,, nA , is a dominant phase 
combination. Therefore we must have i = e n , meaning that A[i] is a pivot column. 

Part (ii) follows from the rank conditions in Lemma I3.4l ii. In particular, rank(<2,j|A[/]) = 
rank(<2/j|A[j, j]) = s+l implies that A[j] € span(A[/], . . . ,A[j — 1]). Since A is in RREF, none of its pivot 
column can be spanned by the preceding columns. Hence A [j] cannot be a pivot column. □ 
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Lemma 1331 identifies outgoing and incoming asymptotic line soli tons respectively with the pivot and 
the non-pivot columns of A. It is then natural to ask if in fact each of the N pivot columns and each of the 
M — N non-pivot columns identifies an outgoing or incoming line soliton, and whether such identification is 
unique. Both of these questions can be answered affirmatively by the following theorem which constitutes 
one of the main results of this work, and is proved in the Appendix. 

Theorem 3.6 (Asymptotic line solitons) Let %fj^{x,y,t) be the tau-function in Eq. (12.11) associated with a 
rank N, irreducible coefficient matrix A with non-negative minors. 

( i) For each pivot index e n there exists a unique asymptotic line soliton as y — > °°, identified by an index 
pair [e n ,j n ] with n = l,...,N and 1 < e„ < j n < M. 

(ii) For each non-pivot index g n there exists a unique asymptotic line soliton as y — > — °°, identified by an 
index pair [i n ,g n ] with n = 1,.. . ,M — N and 1 < i n < g„ < M. 

Thus, the solution ofKPII generated by the coefficient matrix A via Eq. (I2.lt has exactly N + = N asymptotic 
line solitons as y — > °° and N- = M — N asymptotic line solitons as y — > — oo. 

Part (i) of Theorem 13. 61 uniquely identifies the asymptotic line solitons as y — ► oo by the index pairs [e n ,j n ] 
where e„ < j n . The indices e\,...,e^ label the N pivot columns of A, however, the j„'s may correspond 
to either pivot or non-pivot columns, and indeed both cases appear in examples. Moreover, when the pivot 
indices are sorted in increasing order 1 = e\ < £2 < • • • < ?n < M, the indices 71, . . . Jn in general are not 
sorted in any specific order. For example, the line solitons as y — > 00 generated by the matrix A in Eq. (14.51) 
of section|4]have j\ < 73 < }%. In fact, the indices j\ , . . . , need not necessarily even be distinct. Similarly, 
part (ii) of Theorem 13.61 uniquely identifies the asymptotic line solitons as y — > — 00 by index pairs [i n ,g n ], 
where i n < g n . In this case, the indices gi,... ,gM—N label the M — N non-pivot columns of A, but the i„'s 
may correspond to either pivot or non-pivot columns. Moreover, when the non-pivot indices are sorted in 
increasing order I < gi < ■ ■ ■ < gM-N = M), the indices i\ , . . . , im-n are not in general sorted, and need not 
be distinct. Theorem 13.61 yields an important characterization of the solution via the associated coefficient 
matrix A.It provides a concrete method to identify the asymptotic line solitons as y — > ±00, as illustrated 
with the two examples below. Further examples are discussed in section |4] 

Example 3. 7 Consider the tau-function x^m with N = 2 and M = 5 generated by the coefficient matrix 



The pivot columns of A are labeled by the indices {e\ , ei } = { 1 , 3}, and the non-pivot columns by the indices 
{gi,g2,g3} = {2,4,5}. Thus, from Theorem 13.61 we know that there will be = N = 2 asymptotic line 
solitons as y — > 00, identified by the index pairs and [3, j'2] for some j\ > 1 and j% > 3, and that there 

will be Af_ = M — N = 3 asymptotic line solitons as y — > —00, identified by the index pairs [z'i , 2], [12, 4] and 
[h,5], for some ii < 2, ii < 4 and i% < 5. We first determine the asymptotic line solitons as y — > 00 using 
part (i) of Theorem 13. 6l together with the rank conditions in Lemma l3~4l i. Then we find the asymptotic line 
solitons as y — > —00 using part (ii) of Theorem l3.6l and the rank conditions in Lemma l3~4l ii. 

For the first pivot column, e\ = 1, we start with j = 2 and consider the submatrix = (? / j 2 ) ■ Since 
rank(Pi2) = 2>l=Af — 1, from Lemma l3~4l i we conclude that the pair [1 , 2] cannot identify an asymptotic 
line soliton as y — > 00. Incrementing j to j = 3,4,5 and checking the rank of each submatrix P\j we find 




1 
1 




(3.6) 
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Figure 3: Line soliton solutions of KPII: (a) the (3,2)-soliton solution generated by the coefficient matrix A in Exam- 
ple !3.7l with (k\,... ,£5) = (—1,0,4, 1 , |) at t = —32; (b) the inelastic 3-soliton solution generated by the coefficient 
matrix A in Example l3.8l with (ki, . . . ,k^) = (— 1,— A,0, |,1,|) at t = 20 (see text for details). 

that the rank conditions in Lemma l3~4"l i are satisfied when j = 4: = (~, 2 ) = A [5], so rank(Pi 4) = 1 and 
rank(Pi 4|A[l]) = rank(Pi,4)|A[4]) = 2 (The condition rank (Pi ^^[1,4]) = 2 is trivial here, since any three 
columns are linearly dependent.) Thus, the first asymptotic line soliton as y —> 00 is identified by the index 
pair [1,4]. For the second pivot, e 2 = 3, proceeding in a similar manner we find that j = 4 does not satisfy 
the rank conditions, since P 3i 4 has rank 2. But j = 5 satisfies Lemma l3~4l i. since P3 5 = ( j ^ 2 ) , which 
yields rank(P 3 , 5 ) = 1 and rank(P 3 , 5 |A[3]) = rank(P 3 , 5 )|A[5]) = 2. (Again, rank(P 3 , 5 |A[3,5]) = 2 is trivially 
satisfied here.) So the asymptotic line solitons as y — > °° are given by the index pairs [1,4] and [3,5], and the 
associated phase transition diagram (cf. Corollary 12. 61) is given by 

"1,3 >K 1,5 ^4,5 • 

We now consider the asymptotics for y — > —00. Starting with the non-pivot column g\ = 2, the only 
column to its left is i = l. We have 61,2 = 0, andrank(ei i2 |A[l]) =rank(gi i2 |A[2]) =rank(gi j2 |A[l,2]) = 1. 
Consequently, the pair [1,2] identifies an asymptotic line soliton as y — > —00. For g 2 = 4 we consider i = 
1,2,3 and find that the rank conditions in Lemma l3~4l ii are satisfied only for i = 2: in this case, j2 2 ,4 = ( ? ) = 
A[3], so rank(2 2j4 ) = 1 =iV-l and rank(<2 2 , 4 |A[2]) = rank(2 2i4 |A[4]) = 2, while rank(2 2i4 |A[2,4]) = 2 is 
trivially satisfied. Hence [2,4] is the unique asymptotic line soliton as y — > —00 associated to the non-pivot 
column g2 =4. In a similar way we can uniquely identify the last asymptotic line soliton as y — > —00 as 
given by the indices [3,5]. The phase transition diagram for y — > —00 is thus given by 

1— >2 2^4 3^5 
"1,3 ^2,3 ^"3,4 > K$,5 ■ 

To summarize, there are = 2 outgoing line solitons, each associated with one of the pivot columns e\ = \ 
and e 2 = 3, given by the index pairs [1,4] and [3,5], and there are = 3 incoming line solitons, each 
associated with one of the non-pivot columns gi = 2, g 2 = 4 and g 3 = 5, given by the index pairs [1,2], [2,4] 
and [3,5]. A snapshot of the solution at t = —32 is shown in Fig. [3^. 
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Example 3.8 Consider the tau-f unction with N = 3 and M = 6 generated by the coefficient matrix in RREF 

'1110 
,1 | 1 -1 | . (3.7) 
v 1 2 

Again, we first determine the asymptotic line solitons as y — > oo; then we find the asymptotic line solitons as 

y — > — °°. 

The pivot columns of A are labeled by the indices e\ = 1, e% = 4 and e 3 = 5. Thus, we know that the 
asymptotic line solitons as y — > oo will be given by the index pairs [1, ji], [4,72] and [5,73] for some j\, ... ,73. 
Starting with the first pivot, e\ = 1, we take j = 2,3,... and check the rank of the submatrix P,-j in each 

case. When 7 = 2 we have Pi 2 = fo 1 -1) , so rankfPi 2) = 3 > N — 1. So, by Lemma l3~4l i. the index 

Vo 1 2 / ' 

pair [1,2] does not correspond to an asymptotic line soli ton as y — > 00. (In fact, using Lemma l3~71 it can be 
verified that 63,5,6 is the only dominant phase combination along the line 61 = 62 as y — > 00.) We then take 

7 = 3: in this case we have Pi 3 = (00 -lj , with rank(P 1 , 3 ) = 2 =: r andrank(Pi 3 |A[1]) = rank(Pi, 3 |A[3]) = 

rank(Pi 3 |A[1,3]) = 3 = r+l. So the rank conditions in Lemma l3~4l i are satisfied. Therefore the index 
pair [1,3] corresponds to an asymptotic line soliton as y — > 00. Moreover, by considering 7 = 4,5,6 one can 
easily check that the rank conditions are no longer satisfied. Thus [1,3] is the unique asymptotic line soliton 
associated with the pivot index e\ = 1 as y — > 00, in agreement with Theorem 13 .61 Let us now consider the 
second pivot column, £2 = 4. In this case we find that the rank conditions are only satisfied when 7 = 5, since 
P 4 5 = (000 -1) , with rank(P 4 , 5 ) = 2 =: r and rank(P 4 , 5 |A[4]) = rank(P 4 , 5 |A[5]) = rank(P 4 , 5 |A[4,5]) = 3 = 

r + 1 . Therefore, the index pair [4, 5] corresponds to an asymptotic line soliton as y — > 00. Finally, since 
e 3 = 5 and since we know from Theorem 13 .61 that 7 > e 3 , we immediately find that the third asymptotic line 
soliton as y — > 00 is given by the index pair [5,6]. From Corollary 12.61 the phase transition diagram as y — > 00 
is given by 

"1,4,5 ^1,4,6 ^1,5,6 ► "3,5,6 • 

The non-pivot columns of the coefficient matrix A are labeled by the indices g\ = 2 g2 = 3 and g 3 = 
6. For g\ = 2, the only possible value of i < j is i = 1. In this case Q\ 2 = 0, so rank(2i 2) = and 
rank(<2i,2|A[l]) = rank(<2i,2|A[2]) = rank(<2i,2|A[l,2]) = 1. Thus the pair [1,2] identifies an asymptotic line 
soliton as y — > —00. For g2 = 3 we consider i = 2,1: when i = 2, the rank conditions in Lemma l3~4l ii are 
satisfied, leading to the asymptotic line soliton [2,3] as y — ► —00. We can check that the soliton associated 
with the non-pivot column g2 = 3 is unique by considering i = 1 and verifying that the rank conditions 
are not satisfied. Similarly, it is easy to verify that for g 3 = 6 the index pair [4, 6] uniquely identifies the 
asymptotic line soliton as y — > —00. The phase transition diagram as y — ► —00 reads as follows: 

"1,4,5 ► "2,4,5 * "3,4,5 K "3,5,6 • 

Summarizing, there are = 3 asymptotic line solitons as y — > 00, each associated with one of the pivots 
e\ = 1, ei = 4 and e 3 = 5, and indentified by the index pairs [1,3], [4,5] and [5,6], and there are = 3 
asymptotic line solitons as y —> — oo, each associated with one of the non-pivot columns gi = 2, gj = 3 and 
g3 = 6 and identified by the index pairs [1,2], [2,3] and [4,6]. A snapshot of the solution at t = —20 is shown 
in Fig.|3j5. 

Examples 13.71 and 13.81 illustrate the fact that, starting from any given coefficient matrix A in RREF, 
the asymptotic line solitons as y — > ±00 can be identified in an algorithmic way by applying Theorem 13.61 
together with the rank conditions in Lemma l3~4l 
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4 Further examples 



In this section we further illustrate the asymptotic results derived in sections and [3]by discussing a variety 
of solutions of KPII generated by the tau-function dl.3t with different choices of coefficient matrices. 

Ordinary Af-soliton solutions. These are constructed by taking M = 2N and choosing the functions 
{fn}n=i in Eq- CEP as ( e -8-> see Refs - Elini) 

Ux,y jt )=e^"- l +e^", n = l,...,N. (4.1) 

The corresponding coefficient matrix is thus given by 

/I 1 ••• 0\ 
1 1 ••• 

A= : : : : : : : ' 
\0 ••• 11/ 

with N pairs of identical columns at positions {2n — l,2n},n = l,...,N. There are only 2 N non-zero minors 
of A, which are given by A(m\,m2, ■ ■ ■ ,m^) = 1 where, for each n = 1, .. . ,N, either m n = 2n — 1 or m n = 2n. 
The asymptotic analysis presented in the previous section allows one to identify these solutions as a subclass 
of elastic Af-soliton solutions. More precisely, the N solitons are identified by the index pairs [2n — 1 , 2n] 
torn = 1, . .. ,N, where i n = 2n — \ and j„ = 2n label respectively the pivot and non-pivot columns of A. 
Therefore their amplitudes and directions are given by a n = &2„ — &2n-i and c n = ^2«-i +^2«- Moreover, 
the dominant pair of phase combinations for the ra-th soliton as y — > oo is given by 613 2n-i,2n+2,2n+4,...,2N 

and 0i,3,... j 2«-3,2n,2n+2 2N, while the dominant phase combinations for the same soliton as y — ► —00 by 

02,4,...,2»,2ra+i,2n+3,...,2#-i an d 62,4,... ,2re-2,2n-i,2n+l,... ,2iv-i- Apart from the position shift of each soliton, the 
interaction gives rise to a pattern of ,/V intersecting lines in the xy-plane, as shown in Fig. |4^. 



Solutions of KPII which also satisfy the finite Toda lattice hierarchy. Another class of (N -,N+) -soliton 
solutions of KPII is given by the following the choice of functions {f n }^ =l in Eq. (11.101) : 



fn=f 



(*-l) 



1, 



,,N. 



(4.2) 



In addition to generating solutions of KPII, the set of tau-f unctions Tjv,m ior N = 1, . . . ,M also satisfy the 



Pliicker relations for the finite Toda lattice hierarchy \2\. Choosing f(x,y,t) = 
following coefficient matrix: 

/ 1 1 ••• 1 \ 



M 



-1 



=l e "• then yields the 



(4.3) 



Note that A in Eq. (14.31) is not in RREF yet, and coincides with the matrix K in Lemma 12.11 Here the 
pivot columns are labeled by indices 1, . . . ,N; all the (^) minors of A are non-zero, and coincide with the 
Van der Monde determinants (12.4b . This class of solutions was studied in Ref. Q, where it was shown that 
the ,/V asymptotic line solitons as y — > 00 are identified by the index pairs [n,n +M — N] for n = l,...,N, 
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while the M—N asymptotic line solitons as y — > — oo are identified by the index pairs [n,n +N] for n = 
1, . . . ,M — N. These pairings can also be easily verified using Theorem 13.61 The dominant pair of phase 
combinations for the «-th soliton as y — > oo is given by 61 n ,M-N+n+i,...,M an d ®i,...,n-l,M-N+n,...,M, while 
the dominant pair of phase combinations for the n-th soliton as y — > —00 by B n ,.,.fl+n-i an d Q n +i,...fl+n- The 
solution displays phenomena of soliton resonance and web structure (e.g., see Fig. |4Jd). More precisely, 
the interaction of the asymptotic line solitons results in a pattern with (2N- — l)N + interaction vertices, 
(3iV_ — 2)N+ intermediate interaction segments and (N- — l)(N + — l) "holes" in the xy-plane. Each of the 
intermediate interaction segment can be effectively regarded as a line soliton since it satisfies the dispersion 
relation (11.8b - Furthermore, all of the asymptotic and intermediate line solitons interact via a collection of 
fundamental resonances: a. fundamental resonance, also called a Y-junction, is a travelling-wave solution of 
KPII describing an intersection of three line solitons whose wavenumbers k fl and frequencies co a (a = 1 , 2, 3) 
satisfy the three-wave resonance conditions lTP71[T9l 

ki + k 2 = k 3 , coi + (O2 = (03 • (4.4) 

Such a solution is shown in Fig.^. 



Elastic A^-soliton solutions. As mentioned in sections ^ and |3j elastic Af-soliton solutions are those for 
which the sets of incoming and outgoing asymptotic line solitons are the same. In this case we necessarily 
have M = 2N. Ordinary Af-soliton solutions and solutions of KPII which also satisfy the finite Toda lattice 
hierarchy with M = 2N are two special classes of elastic Af-soliton solutions. However, a large variety of 
other elastic Af-soliton solutions do also exist. For example, Fig. @J: shows an elastic 3-soliton solution 
generated by the coefficient matrix: 




1 1 



1 



(4.5) 



1 



In this case the pivot columns are labeled by indices 1, 2 and 3. So, from Lemma 1331 we know that the 
asymptotic line solitons as y — ► 00 will be identified by index pairs [l,ji], [2,72] and [3, 73], while those 
as y — ► —00 by index pairs [z'i,4], [12 , 5] and [z'3,6], for some value of ii , . . . , 13 and Indeed, use 

of the asymptotic techniques developed in section |3] allows one to conclude that both the incoming and 
the outgoing asymptotic line solitons are given by the same index pairs [1,4], [2,6] and [3,5]. The soliton 
interactions in this case are partially resonant, in the sense that the pairwise interaction among solitons 
[1,4] and [2,6] and that among solitons [1,4] and [3,5] are both resonant, but the pairwise interaction among 
solitons [2,6] and [3,5] is non-resonant. Similarly, Fig. @Jl shows an elastic, partially resonant 4-soliton 
solution generated by the coefficient matrix 



f\ -1 1 0-1 -2\ 

12 0-101 2 

1 2 -1 -2 

\0 0000 1 2 3/ 



(4.6) 



In this case the pivot columns are labeled by the indices 1, 2, 4 and 6 and the non-pivot columns by the 
indices 3, 5, 7 and 8. The asymptotic line solitons as y — > ±00 are identified by the index pairs [1,3], 
[2,5], [4,7] and [6,8]. As can be seen from Fig. the pairwise interaction of solitons [1,3] and [2,5], 
solitons [2,5] and [4,7], and [4,7] and [6, 8] are resonant, but all other pairwise interactions (e.g., the pairwise 
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interactions between solitons [1,3] and [4,7], [1,3] and [6,8], [2,5] and [6,8]) are non-resonant. It should be 
clear from these examples that a large variety of elastic Af-soliton solutions with resonant, partially resonant 
and non-resonant interactions is possible. The properties of elastic Af-soliton solutions are studied in detail 
inRefs. lEflfTHl. 

Inelastic A7-soliton solutions. N-soliton solutions that are not elastic are called inelastic. We have already 
seen such solutions in Examples l3.7l and l3~8l (cf. Figs.|3k,b) of section[3] As a further example, Fig.|4^ shows 
an inelastic 2-soliton solution generated by the coefficient matrix 

-Co ::?)■ 

In this case the pivot columns are labeled by indices 1 and 2; the asymptotic line solitons as y — > — oo are 
identified by the index pairs [1,4] and [2,3], while those as y — > oo by the index pairs [1,3] and [2,4]. Notice 
that the outgoing solitons interact resonantly (via two Y-junctions), while the incoming soliton pair interact 
non-resonantly. This is in contrast with an elastic 2-soliton solution, where both incoming and outgoing 
pairs of solitons exhibit the same kind of interaction. Similarly, Fig. |4t" shows inelastic 3-soliton solution 
generated by the coefficient matrix 

/\ 0-1-10 2\ 
A = I 1 2 1 -l] . (4.8) 
\0 1 I / 

Here the pivot columns are labeled by indices 1, 2 and 5; the asymptotic line solitons as y — > oo are identified 
by the index pairs [1,3], [2,5] and [5,6], while those as y — > — oo by the index pairs [1,3], [2,4] and [3,6]. 
Finally, in the generic case one has M ^ 2N, and the numbers of asymptotic line solitons as y — > ±oo are 
different, as in the solutions shown in Figs.|3k and|4j). 

We should point out that one-soliton solutions, ordinary two-soliton solutions and fundamental reso- 
nances have the property that their time evolution is just an overall translation of a fixed spatial pattern. 
The same property does not hold, however, for all the other solutions presented in this work. That is, the 
interaction patterns formed by these line solitons, and the relative positions of the interaction vertices in the 
xy-plane are in general time-dependent. 



5 Conclusions 

In this article we have studied a class of line-soliton solutions of the Kadomtsev-Petviashvili II equation by 
expressing the tau-function as the Wronskian of N linearly independent combinations of M exponentials. 
From the asymptotics of the tau-function as y — > ±oo we showed that each of these solutions of KPII is 
composed of asymptotic line solitons which are defined by the transition between two dominant phase com- 
binations with N—l common phases. Moreover, the number, amplitudes and directions of the asymptotic 
line solitons are invariant in time. We also derived an algorithmic method to identify these asymptotic line 
solitons in a given solution by examining the N x M coefficient matrix A associated with the correspond- 
ing tau-function. In particular, we proved that every N x M, irreducible coefficient matrix A produces an 
(N-,N + ) -soliton solution of KPII in which there are N + = N asymptotic line solitons as y — > oo, labeled 
by the pivot columns of A, and N- = M—N asymptotic line solitons as y — > — oo, labeled by the non-pivot 
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columns of A. Such solutions exhibit a rich variety of time-dependent spatial patterns which include resonant 
soliton interactions and web structure. Finally, we discussed a number of examples of such (iV_,Af + )-soliton 
solutions in order to illustrate the above results. 

It is remarkable that the KPII equation possesses such a rich structure of line-soliton solutions which are 
generated by a simple form of the tau-function. In this work we have primarily focused on the asymptotic 
behavior of the solutions as y — > ±00, but not on their interactions in the xy-plane. A full characterization 
of the interaction patterns of the general (N-,N + ) -soliton solutions is an important open problem, which is 
left for further study. Nonetheless, we believe that our results will provide a key step toward that endeavor. 
Solutions exhibiting phenomena of soliton resonance and web structure have been found for several other 
(2+l)-dimensional integrable systems, and those solutions can also be described by direct algebraic methods 
similar to the ones used here. Therefore we expect that the results presented in this work will also be useful 
to study soli tonic solutions in a variety of other (2+l)-dimensional integrable systems. 
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Appendix 

A. 1 Proof of Theorem EH 

To prove part (i) of Theorem 12.51 it is sufficient to show that, along each line L c , the sign of the inequalities 
among the phase combinations in Definition 12.31 remain unchanged in time as y — > ±00. For this purpose, 
note that the sign of W] ,„ N — m / ^ m > in Eq. d2.6bb is determined by the coefficient of y on the right- 
hand side as y — > ±00 and for finite £ and t, if this coefficient is non-zero. For generic values of the phase 
parameters k\,... ,km this coefficient is indeed non- vanishing, and its sign depends only on the direction c 
of the line L c . Consequently, the dominant phase combinations asymptotically as y — > ±00 are determined 
only by the constant c for finite time. 

Part (ii) of the theorem is proved by showing that the only possible phase transitions are those in which a 
single phase, say 6 m changes to 8 m / between the two dominant phase combinations across adjacent regions, 
and that no other type of transitions can occur. We first prove that single-phase transitions are allowed; then 
we show that no other type of transitions are allowed. In the following, we will assume t to be finite so that 
the dominant phase combinations remain invariant, according to part (i). Suppose that mi! ..., mjv is the domi- 
nant phase combination in a region R asymptotically for large values of |y|. Since R is a proper subset of M 3 , 
it must have a boundary, across which a transition will take place from mil ... jfnjv to some other dominant 
phase combination. S ince Qmi,...,mfl is dominant, A(m\, . . . ^ according to Definition 123] Therefore, 
the columns A[m\], . . . ,A[m#] of the coefficient matrix form a basis of M. N , and for all j ^ {mi, . . . ,m^} we 
have that A [7] is in the span of A[mi], . . . ,A[m#]. Thus there exists at least one column A[m s ] such that the 
coefficient of A [m s ] in the expansion of A [j] is non-zero. We then have A (m\ , ... , m s - 1 , j,m s +\ m^) ^ 0, 
implying that the phase combination Q mu ... Jns - 1 .j,m s+l ,...,m N is actually present in the tau-function. Then, for 
any j ^ {mi ,m^} it is possible to have a single-phase transition from R to the adjacent region R' across 
the line Q„ h = Qj, since the sign of 0„ Js — 0; changes across this line, implying that Q mu ...jn ! _ 1 .j. ms+1 ...., mN is 
larger than 8 m , „ N in R'. 
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We next prove that no other type of transitions can occur apart from single-phase transitions; we do 
so by reductio ad absurdum. Suppose that at least two phases 6 m ,,9 m2 from the dominant phase combi- 
nation m [,..., /)Ja , in a region R are replaced with phases 0^,0,^ during the transition from R to an adja- 
cent region R' . This transition occurs along the common boundary of R and R', which is given by line 
L : (0 mi + 0,„ 2 ) — (8 m ' + Q m > 2 ) = 0, Thus, along L, the differences m , — m / and m2 — m ^ (or, equivalently, 
the differences m , — m ^ and m2 — m / ) must have opposite signs or be both zero. 

If both differences are zero along L, the lines m , = Q m > and 0„, 2 = m ^ (or, equivalently, the lines Q m = 
Q m ' 2 and 8,„' = m ,) must both coincide with the line L in the xy-plane. This is possible only at a given instant 
of time and if the directions of the two lines are the same, i.e., if k mi +& m ' ( = k mi + k m ' 2 (or, equivalently, 
Kn +k m * 2 = 1^' +k mi ). So for generic values of the phase parameters, or for generic values of time, this 
exceptional case can be excluded. Hence, we assume that m , — m / and 0,„ 2 — 8 m > are of opposite signs. 
Note however that 8„„ - 8 m < = mi , ...,„,„ - Q m ' um2l ... >mN and m , - 8 m < = mi) ..., mA , - Q mum > 2>m3 ..., mir Moreover, 
both of these phase differences must be positive in the interior of R if the minors A(m' 1 ,nt2, . . . ,m^) and 
A(m\,m' 2 ,mT, . . . ,m^) are non-zero, since m ,,..., mjV is the dominant phase in R. Hence, we must conclude 
that mi — 8 m < and m2 — 8 m / cannot have opposite signs unless one or both of the phase combinations 
m / mj mjv and mi , m2 , m3 ..., mjv is absent from the tau-function. This requires that either A (m\,m2, ■ ■ ■ ,m^) or 
A(m\,m' 2 ,mT, . . . ,m^) must be zero. A similar argument applied to the the phase differences mi — Q m ' 2 and 
0m 2 ~^m' l leads to the conclusion that one or both of the minors A (m' 2 ,m2, ■ ■ ■ ,m^) and A (mi, m\, 1113 . . .,m^) 
must vanish. However, from the Plucker relations among the N x N minors of A we have 

A{m\,m 2 . . ■ ,m N )A(m[,m 2 , . . . ,m N ) = 

A(mi,m 2 ,mj, . . . ,m^)A(m\,m2, ■ ■ ■ ,m^) — A (mi, m\ . . . ,m^)A(m' 2 ,m2-, ■ ■ ■ ,m^) . (A.l) 

Then it follows that either A(m\,. . . = or A(m' v m 2 ,m3, . . . ,m^) = 0. But this is impossible since by 
assumption both minors on the left-hand-side are associated with dominant phase combinations. Thus, they 
are both non-zero. Hence we have a reached a contradiction which implies that as y — > ±00, phase transitions 
where more than one phase changes simultaneously across adjacent dominant phase regions, are impossible. 

A.2 Proof of Theorem El 

First we need to establish the following Lemma that will be useful in proving the theorem. 

Lemma A.l IfPn is the submatrix defined in Eq. (I3.2t and e n labels the n-th pivot column of an irreducible 
coefficient matrix A, then N — 1 < rank(/ 5 ( , n£ , )i+ i) <N,Vn= 1, . . . ,N. 

Proof. Recall that the pivot indices are ordered as 1 = e\ < ei < . . . < e# < M for an irreducible matrix 
A. Then it follows from Definition I2.2l ii that, corresponding to each pivot column A[e n ] of an irreducible 
matrix A, there exists at least one non-pivot column A [/'*], with j* > e n , that has a non-zero entry in its «-th 
row. Hence we have A(e\,. . . ,e„_i, j*,e n+ \, . . . ,e^) ^ 0. This implies that the matrix A[l, . . . ,e„ — l,e„ + 
1, . . . ,M] = (P en e„+\ |A[e„ + 1]) which contains the columns A[e{\,. . . ,A[e„-i],A\j*],A[e n+ i],. . . ,A[e N ], has 
rank N. Thus, the rank of P e „e„+i is at least N — 1, and this yields the desired result. □ 

We are now ready to prove Theorem 13.61 We prove part (i) here; the proof of part (ii) follows similar 
steps. The proof is divided in two parts. First we show that for each pivot index e n ,n= l,...,N, there exists 
an index j„ > e n with the necessary and sufficient properties for [e n ,j n ] to identify an asymptotic line soliton 
as y — ► 00; then we prove that such a j n is unique. 
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Existence. The proof is constructive. For each pivot index e n , and for any j > e n , we consider the 
rank of the matrix P ert j = A[l,2, ...,e n — l,j + l,... ,M] starting from j = e n + 1. When j = e n + 1 we 
have Pg n j = P e „, e „+i, and therefore Af — 1 < rank(P ejijeji+ i) < N from Lemma |A~T1 If rank(p, n . en+ i) = N, 
then Lemma l3~4l i implies that the pair [e n ,e n + 1] does not identify an asymptotic line soliton as y — > oo. 
In this case, we increment the value of j successively from e n + 1, until 1 rank(P en j) decreases from N 
to N — \. Suppose j = is the smallest index such that rank(P eji j t ) = N — 1 and rank(p, n j t |A[/*]) = N. 
We next check the rank of rank(p, M s \A[e n ]). Since rank(P en j, ) = N— 1, two cases are possible: either 
(a) rank(P en j t \A[e n ]) = N or (b) rank(P £ , n j t |A[e„]) =N— 1. We discuss these two cases separately. 

(a) Suppose that rank(P enJ jA[e„]) = Af. By construction we have rank(p, n jjA [/'*]) = N, and since 
N = rank(A) one also has rank(P, n ; ; \A[e n , j*]) = N. In this case we set = j n . It follows from Lemma l3~4l 
that the pair [e n ,j n ] satisfies the necessary rank conditions to identify an asymptotic line soliton as y — > oo. 
Next we show that these rank conditions are also sufficient in order to determine a pair of dominant phase 
combinations in the tau function corresponding to the single-phase transition e„ — * j n . Since rank(P en j n ) = 
Af — 1 , it is possible to choose Af — 1 linearly independent columns A [p\], ■ ■ ■ ,A\p^-\] from the matrix P en j„ so 
that for all choices of linearly independent columns A[l\] , . . . ,A[l^-i] £ Pe„j n one has 2 Q Pu ..., PN _ 1 > 0/ lv . 

as y —>■ oo along the transition line L en j n . Furthermore, since rank(P en j n \A[e n ]) = rank(P en j n |A[j n ]) = A^, the 
minors A(e n ,pi,...,p r ) and A(j n ,p u . . . ,p r ) are both non-zero, and thus Q en ,p u -,PN-\ an dQj„,Pu-,PN-i form 
a dominant pair of phase combinations as y — > oo along the direction of L e „j n . 

(b) Suppose that rank(p, n , ; ; |A[e„]) = N— 1. Since rank(P e)i J J = A^ — 1 by construction, this means 
that A[e n ] € span (P e „j,)- However, since A[e n ] is a pivot column, it cannot be spanned by its preceding 
columns A[l],... ,A[e n — 1]. Hence the spanning set of A[e n ] from P e „j, must contain at least one col- 
umn from A[j* + 1], . . . ,A[M]. In this case we continue incrementing the value of j starting from until 
the pivot column A[e n ] is no longer in the span of the columns of the resulting submatrix P e „j- Let j n 
be the smallest index such that A[e n ] is spanned by the columns of Pe„jjA[j n ] but not by those of P e „j n - 
Then, by construction we have rank(P en j n ) =: r < N — 1, and rank(Pg n j n |A[e„]) = rank(P eji ; |A[j„]) = 
rank(P en j n |A[e„, j n ]) = r + 1. The rank conditions in Lemma l3~4l i are once again satisfied for the index 
pair \e n , j n ] thus found. The sufficiency of these conditions can then be established by following similar steps 
as in case (a). Namely, it is possible to choose a set of linearly independent vectors A[h], . . . ,A[l r ] £ P e „,j„ 
and extend it to a basis {A[e„],A[/i], . . . , A [/,•], A [mi], . . . ,A[m. s ]} of M. N , where A[mi], . . . ,A[m s ] € Q e „j„ and 
r+s = N— 1. We then have A(e„,h, ... ,l r ,mi,... ,m s ) / 0, which also implies A(j n ,l\ , . . . ,l r ,ni\, . . . ,m s ) ^ 
since A[e n ] G span (P en) yJ A [/'„]). As in case (a), we can now maximize the phase combinations over all 
such sets {/i, . . . ,/ r mi, . . . ,m s }, and find a set of indices {pi, . . . ,p r ,q\, . . . ,q s } such that $e m p\ r ..,p T> q\ q s 

and Qj ntPl p r ,qi,...,q s form a dominant pair of phase combinations as y — > oo along the direction of L en j n . 

Summarizing, we have shown that for each pivot index e n , n = 1,2, .. . ,N, there exists at least one asymp- 
totic line soliton [e n ,j n ] with j n > e n as y — > oo. Next we prove uniqueness. 

'Note that a value of j such that rank(P e(l j) =N—1 always exists, since for j = M we have P ent M = A[l, . .. ,e„— 1] whose rank 
is n — 1, since A is in RREF. 

2 The existence of such a set is guaranteed because part (i) of the dominant phase condition 13.11 implies that, as y — » °° in the 
[e n ,jn] direction, the phases corresponding to the index set P €n i are ordered as 6i > 02 > • •• > 9 e _i and Qj n + \ < Qj„+2 <■•■ < &m- 
Then, since rank(P en / ) = N — 1, it is possible to select the top N—l phases from the above two lists so that the corresponding 
columns are linearly independent. 

3 Note that this is possible only for n < N, because when n = N the submatrix P eN ,j for any j > contains the pivot columns 
A[ei], . . . ,A[e/v_i], Hence, rank(P f , A ,. J ) = N — 1 and rank(P ew ;|A[ejv]) = N. Consequently, n = N always belongs to case (a) above 
and not to case (b). 
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Uniqueness. Suppose that [e„,j„] and [e n ,j' n ] are two asymptotic line solitons identified by the same 
pivot index e n as y —> oo. Without loss of generality, assume that j' n > j n , and consider the line soliton 
[e n ,f n ]. Lemma l3~4l i implies that rankfR ,/ \A[ f n ]) = Tank(P en j n \A[e n ,f n ]). Hence the pivot columnA[e„] is 
spanned by the columns of the submatrix (P e f \A[f n ]). But by assumption we have (P e „.j'„\A[j' n ]) C P £n j n , 
since j' n > j„. Hence A[e n ] is also spanned by the columns of P en j n . This however implies that rank(P en j n ) = 
rank(Pg n ; \A[e n ]), which contradicts the necessary rank conditions in Lemma l3~4l i for [e n ,j n ] to identify an 
asymptotic line soliton as y — > oo. Therefore we must have j n = j n >. Thus, it is not possible to have two 
distinct asymptotic line solitons as y — > oo associated with the same pivot index e„. Part (i) of Theorem l3.6l 
is now proved. 

A.3 Equivalence classes and duality of solutions 

In this appendix, we investigate the relationship between two classes of KPII multi-soliton solutions with 
complementary sets of asymptotic line solitons. Note that the KPII equation dl.lt is invariant under the 
inversion symmetry (x,y,t) — > (— x, —y, — /). As aresult, if u(x,y,t) is an (M — N,N)-soliton solution of KPII 
with M — N incoming and N outgoing line solitons, then u(—x, —y, —t) is a (N,M — Af)-soliton solution of 
KPII where the numbers of incoming and outgoing line solitons are reversed. It follows from Theorem 13 .61 
that the solution u{— x, — y, — t) must correspond to some tau-function Tm-w,m(^,3 ; ,0 associated with an 
M — N x M coefficient matrix whose pivot and non-pivot columns uniquely identify the asymptotic line 
solitons of u(— x, —y, —t). Before proceeding further, we introduce the notion of an equivalence class which 
plays an important role in subsequent discussions. Let denote the set of all phase combinations Q mi! ... i m N 
which appear with nonvanishing coefficients in the tau-function x(x,y,t) of Eq. (12. 21 . 

Definition A. 2 (Equivalence class) Two tau-functions are defined to be in the same equivalence class if (up 
to an overall exponential phase factor) the set & is the same for both. The set of (N- ,N + ) -soliton solutions 
of KPII generated by an equivalence class of tau-functions defines an equivalence class of solutions. 

It is clear from the above definition that tau-functions in a given equivalence class can be viewed as positive- 
definite sums of the same exponential phase combinations but with different sets of coefficients. They are 
parametrized by the same set of phase parameters k\ , . . . ,k.M, but the constants 6 m o in the phase 6 m are differ- 
ent. Moreover, the irreducible coefficient matrices associated with the tau-functions have exactly the same 
sets of vanishing and non- vanishing minors, but the magnitudes of the non- vanishing minors are different for 
different matrices. The asymptotic line solitons of each solution in an equivalence class arise from the same 
i —s- j single phase transition, and are therefore labeled by the same index pair [i, j] . Theorem l3.6l then implies 
that the coefficient matrices associated with the tau-functions in the same equivalence class have identical 
sets of pivot and non-pivot indices which identify respectively, the asymptotic line solitons as y — > oo and as 
y —> — oo. Thus, solutions in the same equivalence class can differ only in the position of each asymptotic 
line solitons and in the location of each interaction vertex. As a result, any (N-,N + ) -soli ton solution of KPII 
can be transformed into any other solution in the same equivalence class by spatio-temporal translations of 
the individual asymptotic line solitons. We refer to the two tau-functions iN,M{x,y,t) and "Zm-nm^^^) as 
dual to each other if the solution u(— x, —y, —t) produced by the function Tn,m(~ x i ~ — an d the solution 
generated by Xm-a?,m(^ ) 3 ; ,0 are in the same equivalence class. Note that Xn.m(—x, —y, —t) is not exactly a 
tau-function according to Eq. J2.2t . but it can be transformed to a dual tau-function TM-iv,M v, t) whose co- 
efficient matrix B can be derived from the coefficient matrix A associated with the tau-function Xf]^i{x,y,t), 
as we show next. 
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Since A is of rank N and in RREF, it can be expressed as A = G] P, where In is the N x N identify 
matrix of pivot columns, G is the N x (M — N) matrix of non-pivot columns, and P denotes the M x M 
permutation matrix of M columns of A. We augment A with M — N additional rows to form the invertible 
MxM matrix 

s -{o ,iy 

where O is the [M — N)xN zero matrix and Im-n is the (M — N) x (M — N) identity matrix. Let A' be the 
(M — N) xM matrix obtained by selecting the last M — N rows of (S~ 1 ) T . The rank of A' is M — N, and the 
following complementarity relation exists between A and A': 

Lemma A. 3 The pivot columns of A' are labeled by exactly the same set of indices which label the non-pivot 
columns of A, and viceversa. Moreover, if A is irreducible A' is also irreducible. 

Proof. From Eq. (IA.2t and the fact that P 1 = P T for a permutation matrix, we obtain 

(S-r=( % j° T )p, (A3) 

y — Lr 1M-N J 

which implies that A' = [— G T ,Im-n]P- It is then clear that the pivot columns of A'P _1 are its last M — N 
columns which correspond to the non-pivot columns of AP -1 = [In, G], and viceversa. The same correspon- 
dence between pivot and non-pivot columns also holds for A and A' because the columns of both matrices 
are permuted by the same matrix P 1 . This proves the first part of the Lemma. 

To establish that A' is irreducible, note first from Definition l2.2l that the permutation of columns preserves 
irreducibility of a matrix. Since A is irreducibile, Definition l2.2l implies that all rows or columns of G and G T 
are non-zero. Therefore the matrix A'P -1 = [—G t ,Im-n]> an d hence A', are both irreducible. □ 

Note that although A' is not in RREF, it be put in RREF by a GL(A r ,IR) transformation. Next, we define the 
matrix B which is also of rank M — N and irreducible like A', and whose columns are obtained from A' as 

B[m] = (-l) m A'[m], m=l,...,M. (A.4) 

Then using Eqs. (IA.3t and (IA.4L the minors of A can be expressed in terms of the complementary minors 
of B via (see e.g., Ref. |7|. p. 21) 

A{h , • • - , In) = (-1)° det(P) B(m { m M _ N ) , (A.5) 

where o = M(M + \)/2+N{N + l)/2, and where the indices m\ <mi< ■ ■ • < ttim-n are the complement 
of 1 < h < h < • • 'In m {1)2, •• • ,M}. Furthermore, B plays the role of a coefficient matrix for the dual 
tau-function as given by the following lemma. 

Lemma A.4 (Duality) Iftff^f(x,y,t) is the tau-function associated with an irreducible N x M coefficient 
matrix A, then the matrix B defined via Eq. (IA.4t generates a tau-function XM-N,M( x ^y,t) that is dual 
to x NM (x,y,t). 

Proof. Without loss of generality we choose the tau-function lN,M{x,y,t) associated with the given equiv- 
alence class of solutions such that 6 m o = for all m = 1, ... ,M in Eq. (12.2t . Then, using Eq. (IA.5t we can 
express the tau-function as 

x NM (-x,-y,-t) = (-l)°det(P)e- 6 ' M x' M _ NM (x,y,t) , (A.6a) 
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where 

Tw-ffjf(x>y,t) = I V(h,...,l N )B(m h ...,m M -N)e 6mi '" M - N , (A.6b) 

\<mi<m2<-<mM-N<M 

with V(h, . . . Jn) denoting the Van der Monde determinant as in Eq. (I2.2t and where the sum is now taken 
over the complementary indices mi,... ,m,M-N instead of l\,...,lfj. (The number of terms in the sum re- 
mains the same since (^) = ( M M N ))- It should be clear from Eq. dl.2t that both X^,m{— x, — y, — t) and 
t' m _ n N (x,y,t) in Eq. (IA.6at generate the same solution u(x,y,t) of KPII although x' M _ N M (x,y,t) is not a 
tau-function as given by Eq. (12.2b . Moreover, all the non-zero minors of B have the same sign, which is de- 
termined by the sign of (— l)°det(P) > 0. Thus, by replacing each Van der Monde coefficient V{1\, . . . 
by V (mi, . . . ,rriM-N) inEq. (IA.6bt . l' M _ N M (x,y,t) can be transformed into a new tau-function Xif-ff ) u(x,y,t) 
associated with the irreducible coefficient matrix B. Since both x' M _ N M (x,y,t) and XM-Nja (x,y,t) are sign- 
definite sums of the same exponential phase combinations, they generate solutions that are in the same 
equivalence class. Therefore, the tau-functions x NM (x,y,t) and XM-Nj^( x ,y,t) are dual to each other, thus 
proving the lemma. □ 

By applying Lemma IA.41 it is easy to show that part (i) of Theorem 13.61 implies part (ii) and viceversa. 
For example, by applying part (i) of Theorem 13.61 to the tau-function XM-N,M(x,y,t) in Lemma |A~4*1 one 
can conclude that as y — > °°, XM-N,M( x ,y, generates a solution with exactly M — N line solitons, identified 
by the pivot indices gi,... ,gM—N of the associated coefficient matrix B. A Since XM-N^t{x,y,t) is dual to 
XM,N{x,y,t), the solution generated by XM-Njn(x,y,t) is in the same equivalence class as u{— x, — y, — t). 
Consequently, the asymptotic line solitons of u(—x, —y, —t) as y — ► °°, are labeled by exactly the same 
indices g\, . . . ,gM-N- Then it follows that as y — > — °°, there are M — N asymptotic line solitons of the 
solution u{x,y,t) generated by Xifjti (x,y,t). Furthermore, these line solitons are labeled by the same indices 
gi, . .. ,gM-N which are the non-pivot indices of the coefficient matrix A of the tau-function x^ t M{x,y,t),thus 
proving part (ii) of Theorem 13 .61 Similarly, one could also prove part (i) of the Theorem using part (ii) and 
Lemmas IA.4I 

Another consequence of Lemma lA~4l is that the dominant pairs of phase combinations for the asymp- 
totic line solitons of XM-N,N(x,y,t) as y — > oo are the complement of those for the asymptotic line solitons 
of the dual tau-function %Nji(x,y,t) as y — > — °°. Thus, if the dominant pair of phase combinations for 
Xm-n,m( x ^ { ) as y ^ oo along the line L u is given by 6,-, m2; ..., mM _ w and Q jm ,..., mM _ N , the dominant phase 
combinations for x^ t M{x,y,t) asy — > — oo along L,j are 0,y 2 and Qj,i 2 ,...,i N , where the index set {fa, ■ . . ,/jv} 
is the complement of {i,j,m2, ■ ■ ■ ,mM-N} in {L • • • ,M}. 

A particularly interesting subclass of (A^_,A^ + )-soliton solutions is obtained by requiring the solutions 
u(x,y,t) and u{— x, — y, — t) to be in the same equivalence class which is generated by "self-dual" tau- 
functions. These are the elastic A^-soliton solutions of KPII, for which the amplitudes and directions of 
the Af incoming line solitons coincide with those of the ,/V outgoing line solitons, as mentioned in section ^ 
Thus, the set of incoming line solitons and the set of outgoing line solitons can both be labeled by the same 
index pairs { [i n , j n ] y . Clearly, in this case we have = = Af and M = 2N. The detailed properties of 
the elastic A^-soliton solution are studied in Refs. l3l [T3l . Here we only mention one result which is a direct 
consequence of Theorem 13 .61 and the above discussions: 

Corollary A.5 A necessary condition for a set of index pairs {[in,jn]}n=i to describe an elastic N-soliton 
solution is that the indices i\, . . . , ipj and ji, . . . , form a disjoint partition of the integers 1 , . . . , 2N. 

4 Since the ordering of the pivot and non-pivot columns of B is reversed with respect to that of A, if with i < j labels an 
asymptotic line soliton generated by Tm-,/v,m as y — > °°, then j is the pivot index, not i. 
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Proof. From part (i) of Theorem l3.6l the indices i\ , . . . , % for the N asymptotic line solitons as y — > oo label 
the pivot columns of A, and from part (ii) of Theorem 13 .61 the indices ,]n for the N asymptotic line 
solitons as y — > — oo label the non-pivot columns of A. In order for the N asymptotic line solitons as y — > — oo 
to be the same as those as y — > oo, however, the index pairs [/„,./„] must obviously be the same as y — > ±00 
for all n = 1, . . . ,A r . But the sets of pivot and non-pivot indices of any matrix are of course disjoint; hence 
the desired result. □ 

Note however that the condition in Corollary IA.5I is necessary but not sufficient to describe an elastic N- 
soliton solution. It is indeed possible to have A-soliton solutions where the index pairs labeling the asymp- 
totic line solitons as y — > ±00 form two different disjoint partition of integers {1,2,..., 2N}. Such A-soliton 
solutions are not elastic. See, for example the 2-soliton solution in Fig.|4^. 
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Figure 4: Line-soliton solutions of KPII: (a) an ordinary 3-soliton solution with {k\, ... ,k&) = (— 3, —2,0, 1, |,2) at 
t = 4; (b) a fully resonant (3,2)-soliton solution with (ki,...,ks) = (— 1,0, 4, 1, |) at t = —32; (c) an elastic, partially 
resonant 3-soliton solution with A given by Eq. J4.5i and (ki,... ,k^) = (— 1,-1,0,4, §, V) at t = —20; (d) an elastic, 
partially resonant 4-soliton solution with A given by Eq. J4.6I > and (ki,...,k$) = (— 2, — |, — 1, — j,0, j, 1, |) at f = 20; 
(e) an inelastic 2-soliton solution with A given by Eq. (14. 7> and (k\ , . . . ,fct) = (— 1, — ^, 5,2) at f = 16; (f) an inelastic 
3-soliton solution with A given by Eq. ( 14. 8> and (fci, . . . = (—1,-4,0, 4,1, |) at t = 32. 
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